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Within a cosmological framework, we provide a Hamiltonian analysis of 
the Mixmaster Universe dynamics on the base of a standard Arnowitt-Deser- 
Misner approach, showing how the chaotic behavior characterizing the evo- 
lution of the system near the cosmological singularity can be obtained as the 
semiclassical limit of the canonical quantization of the model in the same 
dynamical representation. 

The relation between this intrinsic chaotic behavior and the indeterministic 
quantum dynamics is inferred through the coincidence between the micro- 
canonical probability distribution and the semiclassical quantum one. 



1 Introduction 



The simplest and most interesting generalization of the Friedmann-Lemaitre- Robertson- 
Walker (FLRW) cosmology is the Bianchi IX model, whose geometry has the homogene- 
ity constraint but the dynamics makes allowance of anisotropic evolution of different 
(linearly independent) spatial directions. 

Belinski, Kalatnikov and Lifshitz (BKL) at the end of Sixties derived the oscilla- 
tory regime characterizing the asymptotic evolution near the cosmological singularity, 
describing via a discrete map the resulting chaos. 

This Mixmaster |2] dynamics, whose geometry is invariant under the SO (3) group, al- 
lows the line element to be decomposed as a FLRW model plus a gravitational waves 
packet ([3], jl]) sufficiently far from the singularity. 

The various approaches in terms of continuous variables (i.e. construction of an in- 
variant measure for the system |Sl EI and study of the system covariance [Zj-[II]) have 
shown how this chaotic feature is invariant with respect to any choice of the temporal 
gauge and how in this sense it is intrinsic. 

Such deep nature of the Mixmaster deterministic chaos and its very early appearance in 
the Universe evolution lead to believe in the existence of a relation with the quantum 
behavior the system performs during the Planckian era. 

The aim of this paper is to give a precise meaning to this relation by constructing 
the semiclassical limit of a Schoedinger approach to the canonical quantization of the 
Arnowitt-Deser-Misner (ADM) dynamics whose corresponding probability distribution 
coincides with the (deterministic) microcanonical one. 

In Section |21 are outlined the stochastic properties of the Mixmaster dynamics, being 
isomorphic to a billiard on a Lobachevsky plane and is given a stationary statistical 
distribution within the framework of the microcanonical ensemble. 
In Section El we show the existence of a direct correspondence between the classical and 
quantum dynamics outlined by the common form of the continuity equation for the sta- 
tistical distribution and the one for the first order approximation in the semiclassical 
expansion. 

As a remarkable feature, both formalisms (the classical and the quantum one) are con- 
structed in a generic temporal gauge. 



As well known (see ^M,), the dynamics of the Bianchi type IX model, in terms of generic 
Misner-Chitre-like (MCI) variables (r, ^, 6) and corresponding coniugate momenta, is 
summarized by the variational principle 



2 Statistical Mechanics Approach 




(1) 



2 



where 



Uadm = v/i^TT7 , = {e - l) Pi^ + 



(2) 



and [/(r, ^, 6) denotes a potential term which, asymptotically to the singularity, is mod- 
eled by the potential wall 



= Boo {Hi (e, 6)) + Goo {H^ (e, 9)) + (i^n (e, ^)) 
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being Ha {a = I, m, n) the anisotropy parameters |,10,^, such that the motion of the point 
Universe is restricted in the domain Th (see Figure Q). Such a dynamical scheme has 
the relevant feature to let free the choice of the lapse function unless explicitly chosen 
the form of the function F. The bounces against the potential walls and the instability 
of the geodesic flow on the Lobachevsky plane let the dynamics acquire a stochastic 
feature. Moreover in Th the ADM Hamiltonian 




Figure 1: Domain Th, where the anisotropy parameters Ha are simultaneously greater 
than 0. 



becomes (asymptotically) an "energy-like" integral of motion (ITCadm / dT = 0, say 
'Hadm = e = E = const.. 

As shown in ITUl lllj. the Jacobi metric associated with the principle describes a 
point-universe moving over a Lobachevsky plane reduced by the potential walls Q to 
a billiard. 

The Statistics of the Mixmaster stochasticity is reformulated following the lines pre- 
sented in jH], [Hj and ^1] to derive the invariant measure. 

This system is well-described by a microcanonical ensemble, whose Liouville invariant 
measure Woo can be obtained as a solution of the continuity equation 



.dwoo , sin0 dwoo ^sin0 dw^o 



(4) 
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Over the reduced phase space^ 9} (g) S^, Woo exphcitly reads 

I V 0r„®sj 

Using generic MCI variables, the above invariant measure is independent of the choice 
of the temporal gauge, i.e. of the lapse function. 

In order to restrict the invariant measure to the two-dimensional space {C,,0}, we 
stress how the original Hamiltonian equations in the asymptotic limit for which U 
Uoo =^ e = E = const, get in Th the free geodesic motion [Hj 

d$, dr /—^ dO dr sin (p dcj) dV ^ sin 

rf^" rf^V^ -1COS0, — = — -^=^, _ = __-^=^. (6) 

By this system, we easily rewrite the stationary continuity equation in the form 

d4>Woo + ^^?^oo + ^dgw^ = . (7) 

When assuming the independence of the distribution function Woo from then the 
normalization condition leads to the restricted function 

(e, 6) = r w^{^, e, 4>)d<P = 2™oo(e, 0) . (8) 
Jo 

Thus, we get the following continuity equation for g^o 

dgoo , sin0 dg, 



where plays the role of a parameter and the corresponding microcanonical solution on 
the whole configuration space {^, 9} reads 

i V {^,e}^TH 

We conclude observing how the Hamilton equations retain the same form even when 
the potential walls have a finite (non-zero) value. 



3 Semiclassical Limit of the Quantum Dynamics 

What above outlined is a Mixmaster intrinsic feature and not an effect induced by a 
particular class of references: the whole MCI formalism and its consequences has been 
developed in a framework free from the choice of a specific time gauge. 

^Sl denotes the (p-ciide of the momentum space variable. 
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Since chaos appears close enough to the Big Bang, we infer that it has some relations 
with the indeterministic quantum dynamics the model performs in the Planckian era. 
This relation between quantum and deterministic chaos is searched in the sense of a 
semiclassical limit for a canonical quantization of the model. 

The asymptotical principle describes a two dimensional Hamiltonian system, which 
can be quantized by a natural Schroedinger approach 

= -j-nADMV , (11) 

being ip = ?/'(r, ,^,^) the wave function for the point-universe and, implementing 'Hadm 
(see (121)) to an operator^, i.e. 

i ^ L 0^9, 

d d 
Pi Pi = ~^^Q^ ' P'> ^ Po = ~^^d0 ' ^^^'^ 

the equation (fTT|) rewrites explicitly, in the asymptotic limit U Uoo, 

.dip 




d 1 d U, 



1/2 

^, (13) 



where we took an appropriate symmetric normal ordering prescription and we left Uoo 
to stress that the potential cannot be neglected on the entire configuration space {C,, 6}. 
Being Uoo equal to infinity out of Th, i' requires as boundary condition to vanish outside 
the potential walls, say 

ipidTH) = 0. (14) 

The quantum equation (fT^ is equivalent to the Wheeler-DeWitt one for the same Bianchi 
model, once separated the positive and negative frequencies solutions [12], with the 
advantage that now r is a real time variable. Since the potential walls U^o are time 
independent, a solution of this equation has the form 



V;(r,e,e) = Ecne-^™^V(e,^) (15) 

n=l 



where c„ are constant coefficients and we assumed a discrete "energy" spectrum because 
the quantum point-universe is restricted in the finite region Th', the position (|TH)l in (|T!?|l 



^The only non vanishing canonical commutation relations are 



ih . 
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leads to the eigenvalue problem 



d 



d 



d 



d 
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(16) 



In what follows we search the semiclassical solution of this equation regarding the eigen- 
value -Eoon as a finite constant (i.e. we consider the potential walls as finite) and only 
at the end of the procedure we will take the limit for Uoo- 

We infer that, in the semiclassical limit when — and the "occupation number" n 
tends to infinity (but nh approaches a finite value), the wave function ip^ approaches a 
function Lp as 

(17) 



Jim ^„(e,^)=^(e,^) 



lim Er^ 



n— >oo 



E. 



'oo 



The expression for Lp is taken as a semiclassical expansion up to the first order, i.e. 



•(^,^)exp <^i- 



where r and S are functions to be determined. 

Substituting (fTHj) in (fT^ and separating the real from the complex part, we get two 
independent equations, i.e. 

2 , / r^^\ 2 



E. 



+ 



/r 



classical term 

d 



d 



where we multiplied both sides by fi^ and, respectively. 



d 



1 d f dS\ 



dU ^^-ide\ del 



^ =0. 



(19) 



(20) 
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In the limit h —>■ the second term of (fT^ is negligible meanwhile the first one reduces 
to the Hamilton- Jacobi equation 



2 . ^g^^2 



+ 



1 



^2-1 \d9 

The solution of 1^21^ can be easily checked to be'^ 

1 



^\ =E^ 



Ie^-J^ 



di + kde) , 



const. . 



(21) 



(22) 



•^The discontinuity of this function on the boundary of Th is due to the model adopted and does not 
affect the probabihty distribution. 
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We observe that (j2H), through the identifications 

dS dS f 

^ = PC , " ^ ^ J ^^^^^ + ^^^^^ ' ^^^^ 

is reduced to the algebraic relation 

{e - l) P,' + ^^Po' = EJ . (24) 

The constraint is nothing more than the asymptotic one Ti-^j^j^ = E"^ = const, and 
can be solved by setting 

dS , dS 



^/p^ cos(/>, — =pe = E^^^"^ -Isincj), (25) 

where (p ^ [0, 27r[ is a momentum-function related to C, and 6 by the dynamics. On the 
other hand, by (j22|l we get 



1 / 

Pf. = (27) 
to verify the compatibility of these expressions with (|^ we use the equations of motion 



^ ctg0 ^ — 1 sin0 = c, c = const. . (28) 
The required compatibility comes from the identification k = EaoC. Since 

hm E^ = [^ wii'S^r'' (29) 

we see by that the solution ip (,^, 6) vanishes, as due in presence of infinite potential 
walls, outside Vu. 

The substitution in (j20|l of the positions (j25|l (when the potential walls become infinite, 
because of chaotic dynamics, behaves like an independent variable and plays here the 
role of a free parameter) leads to the new equation 

dr smd) Or , , 

We emphasize how this equation coincides with Q, provided the identification r = Q^a'-, 
it is just this correspondence between the statistical and the semiclassical quantum 



''indeed, when the "energy" of the system is greater than the finite potential walls, the classical motion 
is a free one. 
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analysis to ensure that the indeterminism of the quantum dynamics for the Bianchi IX 
model approaches the deterministic chaos in the considered limit. 

Any constant function is a solution of (jHUj) . but the normalization condition requires 
r = l/An and therefore we finally get 

r I |2 I |2_ / ^ V {^,9} e Th .^.x 

say the limit for the quantum probability distribution as n — oo and h ^ associated 
to the wave function 

^iT,e,O = ^ii,0) exp|-z^r(r)| = v^exp|z J {p^d^ + ped9 - E^dT)^ (32) 

coincides with the classical statistical distribution on the microcanonical ensemble. 
Though this formalism of correspondence remains valid for all Bianchi models, only the 
types VIII and IX admit a normalizable wave function ip{C,,6), being confined in Th, 
and a continuity equation which has a real statistical meaning. 
Since referred to stationary states ipn{^,9), the considered semiclassical limit has to be 
intended in view of a "macroscopic" one and is not related to the temporal evolution of 
the model E 



References 

[1] Belinski V A, Khalatnikov I M and Lifshitz E M (1970) Adv. Phys., 19, 525. 

[2] Misner C W (1969) Phys. Rev. Lett. ,22, 1071. 

[3] Lukash V N, (1974) Sov. Phys.-JETP, 40, n. 5, 792-799. 

[4] Grischchuk L P, Doroshkevich A G and Yudin V M (1975) Sov. Phys.-JETP, 42, n.6, 
943-949. 

[5] Chernoff D F and Barrow J D (1983) Phys. Rev. Lett, 50, 134. 
[6] Kirillov A A and Montani G (1997) Phys. Rev. D, 56, n. 10, 6225. 

[7] Cornish N J and Levin JJ (1997) Phys. Rev. Lett, 78, 998; (1997) Phys. Rev. D, 55, 
7489. 

[8] Montani G (2000) in: The Chaotic Universe, Gurzadyan V G and Ruffini R (eds), 173, 
World Sci. 

[9] Hotter A E and Letelier P S avaliable 'g r-qc/OOllOOT] 
[10] Imponente G P, Montani G (2001) Phys. Rev. D, 63, 103501. 

[11] Imponente G P, Montani G (2002) to appear on Int. Journ. Mod. Phys. (available 
^gr-qc/0106028). 



8 



[12] Kuchar K (1981) in Quantum Gravity 2, a Second Oxford Symposium, edited by 
C.J.Isham, R.Penrose and D.W.Sciama, Clarendon Press Oxford, 1981, 229. 

[13] Kirillov A A and Montani G (1997) JETP Lett, 66, 7, 475. 



9 



